Fractal Spin Glass Properties of Low Energy Configurations 
in the Frenkel-Kontorova chain 
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We study numerically and analytically the classical one-dimensional Frenkel-Kontorova chain in the 
regime of pinned phase characterized by phonon gap. Our results show the existence of exponentially 
many static equilibrium configurations which are exponentially close to the energy of the ground 
state. The energies of these configurations form a fractal quasi-degenerate band structure which is 
described on the basis of elementary excitations. Contrary to the ground state, the configurations 
inside these bands are disordered. 
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The Frenkel-Kontorova (FK) model O describes a one- 
dimensional chain of atoms/particles with harmonic cou- 
plings placed in a periodic potential. This model was 
introduced more than sixty years ago with the aim to 
study commensurate-incommensurate phase transitions 
jl],Q. However, it was also successfully applied for the de- 
scription of crystal dislocations H , epitaxial monolayers 
on the crystal surface B , ionic conductors and glassy ma- 
terials pHd] and, more recently, to charge-density waves 
H and dry friction pjll| . In addition, the FK model has 
also found its implementation in the investigation of the 
Josephson junction chain [y_2|. The physical properties of 
the FK model are very rich. Moreover, different types of 
interaction between atoms can be effectively reduced to 
the case of FK model and, due to that, this model con- 
tinues to attract the active interest of different research 
groups. 

The ground state of the classical FK model is defined as 
the static, equilibrium configuration of the chain, which 
corresponds to the absolute minimum of the chain po- 
tential energy. More than twenty years ago Aubry dis- 
covered p,0,H that the ground state is unique and is 
characterized by a special regular order of atoms in the 
chain. In fact, the positions of atoms in the chain are de- 
scribed by an area-preserving map which is well known 
in the field of dynamical chaos and which is called the 
Chirikov standard map |y| . The density of particles in 
the FK model determines the rotation number of the in- 
variant curves of the map while the amplitude of the 
periodic potential gives the value of the dimcnsionlcss 
parameter K. For K < K c the KAM curves are smooth 
and the spectrum of long wave phonon excitations in the 
chain is characterized by a linear dispersion law starting 
from zero frequency. On the contrary for K > K c the 
KAM curves are destroyed and replaced by an invariant 
Cantor set which is called cantorus. In this regime the 
phonon spectrum has a gap so that the phonon excita- 
tions are suppressed at low temperature. The effects of 



the cantorus on the dynamical properties of the map were 
discussed in |lO|Jl 7J ] . Later M] , on the example of Ising 
spin model to which the FK model can be approximately 
reduced fll9|l , it has been shown that the ground state 
has some well defined hierarchical structure. The main 
features of this structure are determined by the number 
properties of the dimensionless particle density which is 
given by the ratio of the mean interparticle distance to 
the period of the external field. 

In more recent studies |2Q-E3J the attention was mainly 
concentrated on phonon modes in incommensurate one- 
dimensional chains. Indeed, the phonon modes con- 
tribute to the specific heat of the system and hence they 
are responsible for the heat conduction along the chain 
p2| , p3[ . The propagation and localization of phonon 
modes pQ,Elf have been studied for small vibrations of 
particles around their equilibrium positions in the ground 
state. 

However we would like to stress that for K > K c , be- 
sides the ground state there exist other excited equilib- 
rium configurations, corresponding to local minima of the 
potential, with energies very close to the ground state. To 
our knowledge only few studies were dedicated to excited 
equilibrium configurations, see for example M,fL8U24] ■ I n 
this paper we study the properties of these low energy 
equilibrium configurations and determine the structure of 
their energy spectrum and its dependence on the strength 
of the periodic potential and on the chain length. The 
obtained results show that these configurations are ex- 
ponentially close in energy to the ground state and the 
number of configurations grows exponentially with the 
length of the chain. We also show that these configu- 
rations have interesting fractal properties which we will 
describe in detail. The transition between different con- 
figurations can be understood on the basis of elementary 
excitations which we call "bricks". The numerical and 
analytical study of these elementary excitations allows to 
understand and describe the fractal structure of energy 



bands corresponding to equilibrium configurations. Since 
the excited equilibrium configurations are exponentially 
close to the ground state they will strongly contribute 
to the physical system properties at finite temperature. 
Contrary to the ground state in which atoms form a reg- 
ular structure, in the excited configurations this order 
is partially destroyed and some chaotic feature appears. 
In some sense the existence of an exponential number of 
configurations exponentially close in energy, reminds the 
situation in classical spin glasses pa]. However, contrary 
to the usual spin glass models, the FK model is described 
by a simple Hamiltonian without any disorder. Therefore 
the appearance of exponentially quasi-degenerate config- 
urations in the FK system can be viewed as a dynamical 
spin glass model. Thus the rich variety of properties of 
the FK model can find application in different areas of 
physics. 



I. THE MODEL. 

Let us consider a chain of particles with pairwise elastic 
interactions between nearest neighbors: V(xi,Xi—i) = 
v(Xi — Xi-i) = (xi — Xi-i) 2 /2 . This chain is placed in 
a periodic external field: W(xi) — —Kcos(xi), where 
(without any loss of generality) the period is taken equal 
to 2tt. Therefore the Hamiltonian of the FK model reads: 



H = 



r>2 



[Xi X{—i) 



K cos(xi)]. 



(i) 



Here we have taken the mass of the particles and the 
elastic constant equal to one. At the equilibrium the 
momenta Pi = and in addition 



dH 

dx, 



= -Xi+l 



2xi — Xi-\ + K sin(xi) = 0. 



(2) 



After the introduction of new variables Pi+\ = £j+i — 
Xi, this equation can be written in the form of an area- 
preserving map: 



Pt+i = Pi + Ksm(xi) , x i+ i = Xi +p i+ i , 



(3) 



which is known as the Chirikov standard map |lq] . 

We concentrate our investigation on the case of golden 
mean dimensionless particle density v = (v5— 1)/2. This 
irrational value can be approximated by rational approx- 
imants which form the Fibonacci sequence s n with num- 
ber of particles s and chain length L = 2nr . In this way 
the rational approximants are v n — r n /s n = s n —i/s n , 
with s n — 1,2,3,5,8,13... and the average distance be- 
tween particles is a — 2-Kv n . For the map (H) the param- 
eter v determines the rotation number of the invariant 
Kolmogorov-Arnold-Moser (KAM) curve. At the golden 
mean value of v the KAM curve is analytical and smooth 



destroyed and the transition by the breaking of analyt- 
icity takes place [p|. As a result the invariant curve is 
replaced by a cantorus which forms an invariant fractal 
set in the phase space of the map. For the FK model the 
cantorus corresponds to the ground state with minimal 
energy as it was shown by Aubry ^M [D|jl4| . 

In this paper we restrict ourselves to the case with K > 
K c , when each particle is locked by potential barriers 
of the external periodical field and the whole chain is 
pinned. Stable configurations of the chain correspond to 
minima of the potential energy: 



U({x}) = J2 



(Xi 



Xi-l) 



K cos(xi)}. 



(4) 



The static ground state corresponds to absolute min- 
imum given by Aubry's solution. However, as we will 
see in the next section, there are other local minima of 
the chain potential U({x}) which give equilibrium static 
configurations with energy being very close to the ground 
state. The number of such configurational states N cs 
grows exponentially with the chain length s. 



II. ENERGY SPECTRUM OF EQUILIBRIUM 
CONFIGURATIONS. 

In Fig.[y we present a typical result for the integrated 
number N cs of excited equilibrium configurations ver- 
sus their energy difference, per particle, from the ground 
stateAC/ = (U — Uq)/s where Uq is the energy of the 
ground state. Here we have number of particles s = 89, 
number of wells r — 55 and two values of K, K = 5 
and K — 2. This figure shows that the energy of con- 
figurations form a sequence of narrow energy bands the 
width of which is much smaller than the distance between 
bands, at least in the very vicinity of the ground state. 
At higher energies the band width starts to grow and 
eventually nearest bands almost merge into each other. 
It is interesting to note that the number of states in each 
band is practically independent of K as it is shown by 
dashed lines in Figjlf with the increase of K each band 
is shifted to smaller values of MJ (in logarithmic scale) 
but the number of states in each band is not changed. 

It should be stressed, that even at a moderate value 
of K — 2 the energy spacing between the ground state 
and the first excited configuration band is of the order of 
10~ 13 . If one assumes that in Eq. (j|) a unit of energy is 
~ leV, then this band is already excited at temperature 
T ~ 10~ 9 °K! Hence one may conclude, that the pure 
ground state is practically inaccessible even for a chain 
with less than one hundred atoms. 



for K < K c = 0.971635.. |26|. For K > K c the curve is 




FIG. 1. Integrated number of equilibrium configuration 
states N cs as a function of the energy difference AC/ between 
the energy of configuration U and the ground state energy 
Ug, counted per one particle. Here the number of particles 
is s = 89 and the number of wells r — 55. The are shown 
for K = 5 (6 upper segments) and for K — 2 (5 lower seg- 
ments). Horizontal dashed lines show the border between 
energy bands. All equilibrium configurations are counted. 

The total number of different local minima N cs is enor- 
mous and it grows very rapidly with K. Therefore to 
numerically find all these configurations one needs to use 
special methods. Our approach to this problem is the 
following. First we find the ground state by the gradi- 
ent method developed by Aubry |0,g7|,g8|. Then we find 
the excited equilibrium configurations with the help of 
the Metropolis algorithm (g9|. In this method the sys- 
tem is considered at some properly chosen temperature 
T. At given T we can probe the configurations with 
AC7 < T while the probability to find configurations with 
higher energy is exponentially suppressed. Our imple- 
mentation of the Metropolis algorithm looks as follows. 
We start from a certain configuration {x}j, which cor- 
responds to some local minimum Uj of the chain poten- 
tial energy U({x}j). Then we take randomly one of the 
chain particles and try to move it into one of the neigh- 
bour wells. Next, with the new distribution of particles 
among the wells, we search for a new local minimum 
U. A new configuration with C/j+i = U is accepted if 
exp(— {U — Uj)/T) > £, where £ is a random number ho- 
mogeneously distributed in the interval [0,1], otherwise 
we try a new attempt. Notice that our Metropolis pro- 
cedure uses particles jumps from well to well rather than 
(small) variations of their coordinates. In this way we 
solve the problem of the Peierls-Nabarro barriers flq] and 
obtain a method with good performance. Physically the 
Peierls-Nabarro barriers are not important since we are 
interested in static configurations and not in the transi- 
tion rate between different states. 

In general the space of low-energy configurations can 
be viewed as a set of disconnected islands. Therefore 
there is a danger that, starting near one island, we can re- 
main in its vicinity forever. To avoid this we periodically 
heat/freeze the system. During this process we perform 
the above described iterations with chosen temperature 



T. In this way the system can move from one island to 
another and visit different equilibrium configurations. 

Since the number of equilibrium configurations is expo- 
nentially large (see Fig J^), it is not possible to visit and 
count exactly all of them. However their number can be 
counted approximately with sufficiently good accuracy in 
the following way. In the lowest excited band the number 
of equilibrium configurations is not so large and it can be 
computed exactly. In order to determine the number of 
states in the next band we start from a representative 
sample of configurations, which is in fact a small part of 
their total number in one band. Then with the help of 
Metropolis algorithm described above we determine the 
ratio between the number of configurations inside the 
first and second band. To do this we choose the tem- 
perature value T in such a way that T ~ 10AC/ 2 > AC7i 
where AC/i^ are the excitation energies for the first and 
second band counted from the ground state. From the 
computed ratio we determine with sufficiently good ac- 
curacy the total number of configurations in the second 
band. By iterating this process we determine the total 
number of states in all bands. Moreover, by gradually 
changing the temperature T, this procedure can be eas- 
ily adapted to higher excitation energies when the bands 
begin to merge. This allows to compute the total num- 
ber of equilibrium configurations in the system which, for 
K = 5 is of the order of 10 15 . 

The energy band spectrum for different values of K 
is shown in FigM It clearly shows that the number of 
bands becomes larger for larger K and, in addition, the 
lowest bands approach exponentially the ground state. 
The existence of such bands exponentially close to the 
ground state is related to the specific properties of the 
FK chain in the pinned phase (K > K c ). This phase 
is characterized by a phonon gap A ||,|27J due to which 
any static displacement perturbation 5xi of particle 
«o (corresponding to a zero- frequency "phonon" ) decays 
exponentially along the chain: 6xi ex exp(— A- \i — io\)- In 
fact A is the Lyapunov exponent of the map (|3|) computed 
on the cantorus. This exponential decay of perturbations 
is responsible for the appearance of exponentially narrow 
bands exponentially close to the ground state. 




FIG. 2. Band energy spectrum of equilibrium configura- 
tions versus the chaos parameter K (upper scale) and the 
phonon gap A (lower scale). The bands are marked by filled 
areas corresponding to a given K value. The chain is the 
same, as in Fig.fll: r/s — 55/89.. The dashed curves are given 
by the semi-empirical expression (H). 

In order to describe the band positions as a function 
of system parameters it is convenient to label the bands 
by the index k in order of increasing energy. Then the 
energies of the four lowest bands are well described by a 
simple empirical formula, see Fig. @: 



(AU k ) ~ Cexp(-az/ fc sv7% + A2), 



(5) 



where (AUk) k is the average energy of fc-th (excited) 
band, s is the number of particles in the chain, and the 
nunrerical values of parameters are C rs 1, a fa 0.59, 
(3 fa 0.12. It is rather interesting to note that this simple 
formula describes quite well even the region with small 
values of A < 0.8 (K < 2). At larger K (and longer 
chains) this formula can be replaced by its even more 
simple limiting expression 



(AU k )~Cexp(-aP k s\). 



(6) 



According to Eqs.(||),@, the spacings between the 
bands and the ground state drops exponentially with the 
length of the chain. In Fig.[| we present the dependence 
of the band structure on the number of particles s in the 
chain, for the rational approximants r/s of the golden 
mean v. 



10 



3 



10 



10" 



10 



10 





"S 

t 
13f 

21f 
34 t 
55 


1 ■ ■ 




^ ^--^ ™"~~---~_ " 




*v ""-.^ "— — „ - 


20 




30 
40 


t ^ 
89 t 

144 "-_ " 

t 
233 







50 



100 



150 200 250 



FIG. 3. Dependence of the band energies on the size s of 
the chain for K = 2 (A = 0.7859). Dashed lines are given by 
empirical formula (H), horizontal dotted lines allow to com- 
pare the band positions for different s. 

The results presented in Figs. 0,0 show that the sim- 
ple Semi-empirical formula (ra), shown by dashed lines, 
describes the positions of the bands in the interval of 30 
orders of magnitude. It is interesting to note that bands 
are also ordered in sonre horizontal levels (marked by dot- 
ted lines) which are practically independent of the size of 
the chain. However, the band width and the number of 
states inside the band of the same horizontal level grows 



with the chain size s. In the next section we will see how 
all these features can be understood on the basis of the 
spatial properties of the chain structure. 

Finally, in Fig.0 we show that the energy band struc- 
ture is characterized by fractal properties. Here the third 
excited band for the chain with r/s = 55/89 and K = 4 
is shown with subsequently growing resolution (see mag- 
nification factors on the top of figure boxes). The hierar- 
chical structure of the bands is evident. Such a structure 
becomes deeper and deeper with the increase of the chain 
length s. In the next section we show the origin of this 
structure and develop a simple model to describe it. 
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FIG. 4. Fractal energy band structure for a chain with 
K = 4, s = 89 and r — 55. Four hierarchical levels with 
growing resolution are shown from left to right. The total 
magnification factor in the horizontal direction is shown on 
the top. Here U is the chain energy per particle and U m in 
gives the energy of the most left band in each panel. The 
vertical scale N ca gives the integrated number of equilibrium 
configurations counted from the bottom of the most left band 
in each panel. The vertical magnification is changed in ten 
times from left to right. 



III. SPATIAL STRUCTURE OF EQUILIBRIUM 
CONFIGURATIONS. 

A. Structure of the ground state. 

To analyze the origin of the FK chain hierarchical 
structure let us start with study of its ground state. It 
is very instructive to analyze regularities of particle po- 
sitions inside the wells. In particular, their positions 
modulo the period of the external field are given by the 
broadly discussed hull function [[7| p^JT^ , p7[ | . Its typical 
example is presented in Fig]5|. 
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FIG. 5. The particle position xi mod 2-k versus the particle 
index (multiplied the average distance a — 2irr/s) mod 2tt. 
This hull function is shown for r — 89, s = 144 and K = 2. 

In this plot the bottoms of potential wells correspond 
to Xi Mod 2-7T = and 2-7T. It is easy to see, that a con- 
siderable amount of particles is located very close to the 
bottoms. To render this observation even more signifi- 
cant, the absolute values of deviations from the bottom 
versus the particle number i are plotted in logarithmic 
scale in Fig.[| 
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FIG. 6. Absolute value of the particle deviations from the 
nearest potential well bottom versus particle number along 
the chain. The chain parameters are as in Fig.H. 

We see, that some of particles are at the well bottoms 
with extremely good accuracy! Moreover, the values of 
small deviations are grouped into three well resolved hier- 
archical levels. Separations along the chain for these par- 
ticles are also ordered in some regular way. The two par- 
ticles closest to the bottom |Ax| sw 4.7- 10~ 25 (!) are sep- 
arated by the distances 55 and 89 (the chain is periodic). 
Then 8 particles (including the previous two) whose de- 
viation from the bottoms is |Ax| < 3-10 -6 , are separated 
by distances 13 and 21, see FigH Finally, 34 particles 
whose deviation from bottoms is \Ax\ < 10 _1 are sepa- 
rated by distances 3 and 5. The greater is K, the closer 
these particles are to the bottoms, yet their separations 
along the chain remain the same. By taking a chain with 
longer length (r/s = 144/233, 233/377, 377/610, . . .) one 
can observe subsequent levels of the hierarchical struc- 
ture. 



The fact that some particles are very close to the bot- 
tom of the wells, is very important. Indeed let us assume 
for a moment, that these particles are exactly at the well 
bottoms. This means, that tension forces acting from 
both sides on any such particle, called hereafter as "glue" 
particle, balance each other exactly. Now, let us cut the 
chain at glue particles into fragments, or "bricks". Then 
we can interchange any two fragments of the chain with- 
out changing the chain potential energy. In general, the 
interchanged bricks are different, and we get in this way 
a new configuration with the same potential energy. So 
we may conclude, that we can get a combinatorially large 
number of degenerate configurations in the ground state 
whose number grows exponentially with the length of the 
chain. 

In fact, our glue particles are lying very close to, but 
not exactly, the bottoms of wells. Actually they are 
slightly shifted from the bottoms, and therefore the ten- 
sions / at the ends of different bricks are not the same. 
As a consequence, when we exchange two different bricks, 
each brick's end will be slightly distorted. The distortion 
is proportional to the difference in boundary tensions A/ 
of the nearby bricks. This leads to a local change of the 
chain energy 



AC/ ~ AU b + AU„, 



(7) 



where AC/& ~ (A/) 2 /2 is due to the distortion of 
nearby bricks, and AU g w 7^(Aa;) 2 /2 = (Af) 2 /2K is 
the change of potential energy due to the shift of the 
glue particle between the bricks. We note that, since glue 
particle deviations are exponentially small and hierarchi- 
cally ordered, then the corresponding tension differences 
are also exponentially small and ordered. Therefore the 
energy change caused by bricks permutation depends on 
the level of the hierarchy inside which the permutation 
is done. The lowest level of the hierarchy is built by 
bricks of two types, which consist of two and four par- 
ticles respectively. For brevity let us denote them as 2 
and 4. Then a chain which consist of 8 particles can be 
denoted as g2gi (the letter g stands for a glue particle). 
The tension difference at this level of the hierarchy is 
Af = KAx~ 1CT 1 . 

The next level of hierarchy has bricks 12 = (4g2gA) 
and 20 = {Ag2gAg2gA). The brackets are introduced 
for convenience to denote the form of the brick. The 
tension difference at this level is much smaller: A/ = 
K Ax ~ 10 -5 . Finally, the third level of the hierarchy 
is composed in the similar way: 54 = (2Qgl2g2§) and 
88 = (20gl2g20gl2g20), with the corresponding tension 
difference A/ = K Ax ~ 10 -24 . With increasing the 
particle number the above described process proceeds in 
a similar way. A simple estimate for the tension dif- 
ference valid at any hierarchical level, and for any K, 
can be written as: A/ ~ K exp(— As m , n ), where s m i n is 
the number of particles in the smallest brick at the given 



level of hierarchy and A is the phonon gap which depends 
implicitly on K. Notice that a brick with the addition 
of the glue particle, forms an elementary cell the size of 
which is given by the Fibonacci numbers. For example 
( g 2)=3,(g4)=5, ( g 2g4)=8 etc. 

The composition rules for the brick construction at any 
hierarchical level can be summarized in following way. 
Suppose that a given level of hierarchy is composed by 
two bricks A and B, with the length of A smaller than 
B. Then the bricks A' B' of the next level can be built 



A' = BgAgB, B' = BgAgBgAgBg 



(8) 



In principle, the composition rules just described al- 
low to build the ground state for a chain of any length. 
It is also clear that for long enough chains one does not 
need to search the global minimum of the potential en- 
ergy. Instead, it is sufficient to minimize the energy of 
bricks up to some hierarchical level: any further opti- 
mization goes beyond any reasonable precision. This 
however also means that within the same precision the 
ground state configuration described by Aubry is indis- 
tinguishable from exponentially many disordered excited 
configurations. 



B. Structure of the excited configurations. 

The picture of the ground state described above allows 
also to understand the structure of excited configura- 
tions. However, in this case the structure can be a bit less 
self-evident. To illustrate this, in Fig.[j] we plot particles 
deviations from well bottoms for a configuration from the 
first excited band in the chain shown in FigJ(| The hull 
function for a typical configuration in this band is shown 
in Fig.g(a) . The hull function for a typical configuration 
in the second, third and fourth excited bands (see the 
band structure in Fig.g with s — 144) is shown in FigJ|. 
Contrary to the monotonic hull function of the ground 
state, here the hull function becomes not monotonic and 
one can see the overlap between horizontal plateaus. 
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FIG. 7. Deviations of the particles from potential well bot- 
toms for an excited configuration taken from the first excited 
energy band (see FigH for s = 144). The separation of this 
band from the ground state of Fig.pl is AU± ~ 10~ 20 , which 
implies that the deviations of glue particles from the well bot- 
toms are Ax < 10 -10 (dashed line). 

From Fig.0 we see that for the first two levels of hi- 
erarchy, the deviations of glue particles from the well 
bottom are practically the same as in the ground state 
(see Fig. g). However at the third hierarchical level the 
deviations of two glue particles (below the dashed line) 
become considerably larger than the corresponding ones 
in the ground state (see Fig.pl). 
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FIG. 8. The hull function for a typical equilibrium config- 
uration in the fc-th excited band for K = 2 and r/s = 89/144: 
a)fc = 1, b) k — 2, c) k — 3, d) A; = 4. The energy band 
structure is shown in Fig. H. Compare with the hull function 
of the ground state shown in Fig. U. 

In order to give an unambiguous definition of bricks 
let us remind, that we want to split the chain into bricks, 
which permutations keeps the chain configuration in- 
side the same band. According to Eq.(Q), the energy 
change due to a permutation, produced by the tension 
differences between permuted bricks, can be estimated 
as SU ~ (Ax) 2 . Therefore the deviations of glue par- 
ticles from the bottom between the bricks is restricted 
by the condition Ax < (AUk) 1 ^ 2 , where AUk is the 
band energy counted from the ground state. Taking 
this condition into account, we can write for the con- 
figuration shown in Fig.p] its decomposition into bricks 
as g20gl22 = g20g(20gl2g20gl2g20gl2g20) 1 where the 
expansion of the configuration is shown up to bricks of 
the second level, 12 and 20. As above, by brackets we 
mark the chain fragments which in permutations should 
be considered as a single brick, since their destruction 
results in the energy change exceeding the band width. 

Let us now discuss the properties of the bricks ex- 
pansion on the example of a periodical chain with 
r/s = 89/144 and K = 2 (see FigJ|). The first ex- 
cited band k = 1 has the excitation energy AUi w 
9.16 • 10 -21 and is composed from one configuration 



g20g(20gl2g20gl2g20gl2g20)(here we do not count the 
configurations with a shift along the chain and reflec- 
tion) . It is interesting to note that this configuration has 
a long commensurate fragment (123/144). 

The second excited band k = 2 has energy AU2 < 
1(T 12 . 

It is composed by three configurations: 
g 12g20g20g20gl2g20g20gl2, 12g20g20g20g20gl2g20gl2 
and g20g20g20g20g20gl2gl2gl2. With the configuration 
from the first band they give all possible different com- 
binations of three bricks 12 and five bricks 20 which are 
used in the composition of the ground state. 



The third band k — 3 has the excitation energy AtA 



1CP 8 . This band has too many configurations to be listed 
here. Let us however mention a new phenomenon which 
appears in this band namely a brick "chemical" reaction 
with dissociation of larger elementary bricks of the sec- 
ond hierarchical level: 



20 + 20 -> 12 + 28. 20 + 28 -> 12 + 36. 



(9) 



Note, that a "free radical" 8 coming from dissociation 
20 — ► 12 + 8 is easily captured by other long bricks, so 
that there is a considerable contribution of long commen- 
surate structures. Near the bottom of the band a typ- 
ical configuration is g20gl2g29gl2g20gl2gl2g20, While 
at the top one has g53gl2gl2gl2gl2gl2gl2\ 

The fourth band k = 4 has energy AU4 < 10~ 5 . Here 
we see a dissociation of the bricks from the second hier- 
archical level: 



12 + 20 -> 4 + 28, 12 + 28 



36, 



(10) 



and appearance of elementary bricks 4 from the first 
hierarchical level. Here are some examples of configu- 
rations 

in this band with bricks 4: g20g20g28g20gl2gl2g4g20, 
gl2gl2gl2g28gl2g20g4g36, . . .. 

Further steps in the whole picture are straightforward. 
Now we outline a simple theory which turns our quali- 
tative observations into quantitative predictions for the 
band energy spectrum. 



C. An analytical approach. 

In fact, the construction of bricks is based on the ex- 
istence of an intrinsic small parameter which allows to 
develop a simple rapidly converging perturbation theory. 
Here we outline its main elements. Let us consider the 
FK chain with s particles and fixed ends at xq — and 
x s = 2-nr. Then the largest brick contains n = s — 1 
particles. If the glue particles (i = 0, i = s) are slightly 
shifted from the well bottoms x a ,b <C< 1, then the brick 
energy can be written as 



U^{x ai x b ) 



- U {n) 

— U Q 



An) 



(*), 



fk n 'x a + f™x b + 

+R^4+R^4-T^x a x b , 



(11) 



where Uq = U^ (0, 0) is the unperturbed energy, 
fab and R™1 are tensions and rigidities at the left/right 
ends of the brick, and T^ is the static "transmission" 
factor along the brick with nparticles. If the brick is 



(n) _ An) __ 



ti n > = /(») and itf 



(n) 



R 



(n) _ 



symmetric then /„ 
Ry 1 ' . The key point of the theory is that in the presence 
of a nonzero phonon gap A the transmission factor T^ n > 
is exponentially small: T' n ) ~ exp(— An). Therefore it 
can be very efficiently used as an expansion parameter in 
the calculations of the energy band spectrum. 

Suppose that at some hierarchical level we have two 
elementary bricks A and B, with lengths ha < tib- Ac- 
cording to our rule of brick composition (ph , we can cal- 
culate the energy of the brick A' — BgAgB as: 

U [A '\x a ,x b ) = mm[U B {x a , Xi) + U A (x ll x 2 ) + U B (x 2 ,x b ) 

(12) 

Then, rewriting (O) in the form ( pjj ) we obtain the 
transformation rules for brick parameters R, f and T. In 
the leading order approximation in the small parameter 
T these rules have the form: 



R 



A' 



f 



A' 



= R A 
= f B - 



rpA' _ 



T B (f B -f A ) 

R A +R B +K 

(T B ) 2 



(13) 



(R A + R B + K) 2 ' 

In the same way for B' we obtain 



R B ' 

f B ' 

rpB' 



R B 

cA' 



= r + (a/) 

= (R A +R B +K) 4 



(14) 



where the tension difference (A/)' = f B —f A between 
new bricks A' and B' can be expressed through the brick 
tension difference (A/) = f B — f A as: 



(A/)' = 



B\2 



T A {T B ) 



(R A + R B + K) : 



:(A/)- 



(15) 



To apply these transformation rules one needs to know 
the bricks parameters at the lowest hierarchical level, e.g. 
for bricks 2 and 4. In this case the number of particles 
is small and the expansion (ftW can be performed ana- 
lytically. For the case K = 2 considered above, we get 
T< 2 ) = 0.24, RW = 0.454, T^ = 0.0533, R^ = 0.32, 
and (A/)( 2,4 ) = 0.21. By applying the transforma- 
tion rules to these data we obtain for the bricks of the 
next hierarchical level 12 and 20: T^ 12 ) = 8.86 • 10~ 5 , 
T( 2 °) = 1.796 • 10" 7 , (A/)( 12 - 2 °) = -6.93 • lO" 6 . The 
exact numerical simulation gives T^ 12 ' = 9.82 • 10~ 5 , 
T( 2 °) = 1.799 -10- 7 , (A/)( 12 - 2 °) = -7.16 • 10~ 6 . Starting 
with exact values for bricks 12 and 20 the transformation 
rules give for bricks 54 and 89 results which are correct 
within four digits accuracy! 



Therefore this simple approach can quantitatively ex- 
plain the splitting of the whole spectra into bands. 
Surely, the leading terms in the small parameter T , as 
well as the expansion (O), can be insufficient to repro- 
duce with high accuracy the deep levels of hierarchical 
band structure. To this end one should take into account 
higher order terms. 

The results presented in this section show that the 
number of equilibrium configurations grows very quickly 
with the length of the chain and with the chaos parame- 
ter K. These configurations form bands which are placed 
exponentially close to the ground state. As a result, even 
in a fixed very small vicinity of the ground state, the num- 
ber of configurations grows exponentially with the chain 
length. This fact is illustrated in Fig. || 




FIG. 9. The number of equilibrium configurations N cs 
with excitation energy from the ground state AU < 2 ■ 10~ 8 
as a function of the number of particles s in the chain for 
K = 2. The dashed line shows the fitted exponential depen- 
dence N ca = 1.78exp(0.0554s). 



IV. DISCUSSION AND CONCLUSIONS. 

In this paper we studied the properties of equilibrium 
static configurations in the Frenkel-Kontorova chain in 
the regime of pinned phase characterized by phonon gap. 
This FK model is rather general and finds applications 
not only for commensurate -incommensurate transition 
for atoms placed on a periodic substrate but also in many 
other fields of physics. In addition, near the equilib- 
rium, also the cases with long range interactions between 
atoms can be effectively reduced to the FK model with 
only nearest neighbors interaction. We have shown that 
energies of equilibrium configurations form a hierarchi- 
cal band structure so that exponentially many configu- 
rations become exponentially close to the unique ground 
state. In this respect the FK model has certain simi- 
larities with classical spin glass models which also are 
characterized by existence of exponentially many quasi- 
degenerate states |29] . At the same time in the FK model 
the disorder is absent and the quasi-degenerate configu- 
rations form a fractal sequence of energy bands which in 



a sense can be considered a dynamical spin glass. On 
the basis of extended numerical and analytical investiga- 
tions we determined the low energy excitation inside the 
quasi-degenerate bands which have a form of bricks from 
which the whole chain can be composed. On the basis of 
these results we have shown that while the ground state is 
characterized by regular structure, the low energy excited 
configurations are disordered due to elementary brick dis- 
placements. This means that exponentially close to the 
ground state there are disordered configurations which 
may have rather different physical properties compared 
to the ground state. For example this disorder should 
significantly affects the properties of phonon excitations 
in the chain. The exponential quasi-degeneracy of low 
energy configurations should be also important in the 
case of quantum FK chain when quantum particles can 
tunnel from one configuration to another. These two as- 
pects are related with new interesting physical effects of 
low energy excitations in many-body systems and require 
further investigations p0[ . 
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